Photon-number statistics with Silicon photomultipliers 
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>- ■ We present a description of the operation of a multi-pixel detector in the 



presence of non-negligible dark-count and cross-talk effects. We apply the 



model to devise self-consistent calibration strategies to be performed on the 
very light under investigation. © 2010 Optical Society of America 
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^ ■ 1. Introduction 

In the last decade many efforts have been devoted to the development and characterization 
of photon-number resolving detectors. The motivation stems from the variety of fields of 
applications: as an example, this kind of devices has already been used in particle physics 
experiments [THE], PET systems [3l[7], biomedical research [8HT2] and atmospheric pollution 
measurements p3J[14]. Moreover, there are other interesting fields of application in which 
photon-counting detectors could be exploited, such as LIDAR (light detection and ranging), 
laser scanning microscopy, and fluorescence correlation spectroscopy [15]. In addition, this 
class of detectors play a key role in quantum optics as it offers the possibility to characterize 
nonclassical multiphoton states. In this field they can be considered as a valid alternative to 
optical homodyne tomography [TB] as they can be employed in direct detection schemes |17j . 
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Several concepts and technologies have been proposed that lead to the development of 
detectors such as visible-light photon counters (VLPC) [IE], superconductive transition edge 
sensors (TES) pj], time-multiplexed detectors [20H22] . hybrid photodetectors (HPD) [231EH] 
and Silicon photomultipliers (SiPM) [25] . Irrespective of the concept and design features, 
these detectors may in general be classified in terms of photon detection efficiency, spectral 
response, time development of the signal, dead time, and notable photon number resolv- 
ing capability. As of today, the ideal detector has yet to appear and the optimal choice is 
application specific. 

This paper focuses on SiPMs, detectors featuring unique characteristics that are achieved 
by a rapidly evolving technology. Silicon photomultipliers consist of a high density (by now 
limited to ~ 2000 cells/mm 2 ) matrix of diodes with a common output. Each diode (or cell) 
is operated in a limited Geiger-Mueller (GM) regime, in order to achieve gains at the level of 
10 6 . Quenching mechanisms are implemented to avoid establishing self-sustaining discharges. 
These detectors are sensitive to single photons triggering GM avalanches and can be endowed 
with a dynamic range well above 100 photons/burst. The photon detection efficiency (PDE) 
depends on the sensor design and specification, but it may well exceed 60%. Moreover, SiPM 
are genuine photon-number resolving detectors in that they measure light intensity simply 
by the number of fired diodes. Compactness, robustness, low cost, low operating voltage, and 
power consumption are also added values against traditional photodetectors. On the other 
hand, SiPMs are affected by significant dark count rates (DCR), associated to cells fired 
by thermally generated charge carriers. Moreover, the GM avalanche development is known 
to be associated to the generation of photons [26], which may in turn trigger secondary 
avalanches and result in relevant cross-talk. Whether DCR and cross-talk may be directly 
measured, it is clear that they are folded in the detector response to any signal and need to 
be modelled to properly assess the statistical properties of the light field being investigated. 
This paper reports the experimental validation of two models, on the way to a self-consistent 
characterization of the SiPM response. 

2. Experimental set-up 

The detector response to a weak light field is shown in Fig. [D^a), featuring the sensor output 
signal after a high-bandwidth amplifier with a gain of 50. The different bands in the image, 
obtained in persistency mode, correspond to samples with different numbers of triggered 
cells, i.e. different numbers of detected photons. The photon-number resolving properties are 
also clear in Fig. H^b) that shows the corresponding spectrum as obtained by digitizing the 
amplified current pulse from the detector over a well defined time window, synchronized to 
the light pulse. These data and the results reported in the following were obtained with a 
SiPM produced by Hamamatsu Photonics [15]; more specifically, the main features of the 
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detector are reported in Table HJ We have decided to test a detector with 10 2 cells in order 
to keep DCR and cross-talk at a reasonable level and to maximize PDE. The signal was 
integrated by a charge digitizer V792 produced by Caen [27], with a 12-bit resolution over 
400 pC range; the signal was typically integrated over a 200 ns long time window. The DCR 
and cross-talk can be directly obtained by measuring the frequency of the pulses from the 
sensor, with no illumination, above a well defined threshold. The result of a threshold scan 
of DCR is displayed in Fig. 121 clearly showing that the single cell avalanches stochastically 
triggered by thermally generated carriers can induce, by optical cross-talk, a cascade effect 
resulting in a multiple cell spurious event. Thus optical cross-talk can be measured using: 



where v± is the the dark count rate measured by setting a threshold at a value lower than 
the peak amplitude of the signal corresponding to a single avalanche and is the count 
rate obtained when the threshold is set beyond the single cell signals (see Figure [2]). The 
experimental results reported below were obtained illuminating the sensor with a frequency- 
doubled Nd:YLF mode-locked laser amplified at 500 Hz (High Q Laser Production) that 
provides linearly polarized pulses of ~5.4 ps duration at 523 nm wavelength. Two series of 
measurements were performed, the first one directly on the coherent laser output and the 
second one on the pseudo-thermal light obtained by passing the laser through a rotating 
ground-glass diffuser (D in Fig. [3]) [28] . The light to be measured was delivered to the sensor 
by a multimode optical fiber (1 mm core diameter). The signal digitization was synchronized 
to the laser pulse. 

3. Detector response modelling 

The response of an ideal detector to a light field can be described in a simple way as a 
bernoullian process: 



being n the number of impinging photons over the integration time, m the number of de- 
tected photons and r] < 1 the photon-detection efficiency. Actually, rj is a single parameter 
quantifying detector effects and losses (intentional or accidental) which can be tracked to 
the optical system. As far as SiPMs are concerned, detector effects are due to the quantum 
efficiency, the fill factor and the avalanche triggering probability, namely the probability for 
a charge carrier to develop a Geiger-Mueller quenched discharge [291130] . As a consequence, 
the distribution P m>e \ of the number of detected photons, that is the GM avalanches actually 
corresponding to a detected photon has to be linked to the distribution P njP h of the number 
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of photons in the light under measurement by [3TH33] : 

oo 

Pm,e\ = 51 Pm,n{v)Pn,ph . (3) 
n=m 

It can be demonstrated [241134] that for a combination of classical light states the statistics 
is preserved by the primary detection process. This simple description has to be further 
developed to link P m ,ei to the probability density distribution of the GM avalanches of any 
origin. First we must take into account spurious hits and cross-talk effects, not negligible in 
the detectors being studied. The dark count rate results in a poissonian process which can 
be described as: 



P m ,dc = ™ dc / m! exp(-m dc ) , (4) 
where m dc is the mean number of dark counts during the gate window (or integration time) 

nrif ] _(2) _ _(3) _ — 

ana o mdc — cr m dc — m^. 

As a consequence, the statistics of the recorded pulses may be described as: 

Pm,c\+&c ^ ] Pi,dcPm—i,e\ > (5) 
i=0 

obviously shifting the mean value to m e i +dc = m e i+m dc , with an increased variance and third- 
order central moment a^ cl+dc = a^ cl + a%] dc = a® A + m dc and a^ el+dc = <r® d + a^ dc = 

(3) , — 
0m,el + m dc- 

As a further step, cross-talk effects shall be taken into account. Cross-talk is a genuine 
cascade phenomenon that can be described at first order as [35J 

CkM = | k [ Je^(l-ef- fe . (6) 

being e the (constant) probability that the GM avalanche of a cell triggers a second cell, 
I the number of dark counts and photo-triggered avalanches and k the actual light signal 
amplitude. Within this first-order approximation, the actual sensor response is described by 

k 

Pk, cross ^ ] Cfc,m(^)-/"m,el+dc , (7) 
m=0 

characterized by k cross = (1 + e)m e i +dc , cr^ ross = (1 + e) 2 o"^ el+dc + e(l - e)m m>el+dc and 
4^ross = (1 + e ) 3 ^Sci+dc + 3e(l - e 2 )cT^ cl+dc + e(l - 3e + 2e 2 )m miC i +dc . In the following we 
refer to this analytical model as Model I. We also consider a better refined model (vide infra, 
Model II) offering, in principle, an extended range of application, but paying the price of 
being limited to a numerical rather than analytical solution. Irrespective of the model, the 
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amplification and digitization processes that produces the output x can simply be described 
as a multiplicative parameter 7: 



x,out 



7-P7AV 



conveniently scaling the momenta as x out = ^k cross , a^ out = 7 2 4^ ross and a^ out = 7 3 o"£Lss- 

3. A. Model I: an analytical evaluation of the second and third order momenta 

This approach extends the method presented in [23] to detectors with a significant dark 
count rate and first order cross-talk effects. Experimentally, it is based on the detection of a 
light field performed by varying, in a controlled way, the optical losses, i.e. 77, with detector 
parameters presumed to be constant throughout the rj scan. The second-order momentum 
of the recorded pulse distribution P X)(m t can be used to evaluate the Fano factor: 

rr (2) (T {2) 
j-1 u x,aut k, cross 

7" 



cross 



■^out k 

(2) 

7(1 + e) = +7- 



m cl+dc 1 + e 

Qci+dc_ 1 + 3e 

= = Xout + 1— — , (9) 

m el+dc 1 + e 

where Q e i + d c = ^mei+dc/^ei+dc — 1 is the Mandel factor of the primary charges. Note that, 
due to dark-counts, the coefficient of x out in Eq. cannot be written as Q p h/n [23ll36j . 
that is, the coefficient Qei+dc/^ei+dc does not only depend on the light field to be measured. 
Similarly we can calculate a sort of symmetry parameter 

rr {3) cr (3) 

rr u x,out 2 k, cross 

*^a;,out 7 ~ z f z 

•^out rCcross 

2/i , \2 (T rn,el+dc . Q 2 l-\ \ ^m.el+dc . 2 ~~ ^6 + 2e 2 ) 



7 (1 + e) ^ V 37 e(l - e)^ V 7 



m e i+dc "^ci+dc 1 + e 

Qs,cl+dc ~~ 3Qcl+dc_2 1 + 3e Qcl+ dc _ ^ /,m 

= =2 x out + 7^ — = ^out + 7 -7- — , 10 

™el+dc 1 + 6 "icl+dc 1 + £ 

where Q Sj ei+dc = ^mei+dcA^ei+dc ~ 1- Note that in the presence of dark-counts both coef- 
ficients (Q SiG i + dc — 3<5ci+dc)/ ; "Tei+dc an d Qci+dc/^ei+dc are no-more independent of the light 
under measurement [21]. We will see how this modifies the results in the following examples 
involving classical light states for which the statistics of detected photons is the same as that 
of photons. 

A validation and a comparison of the proposed model was performed by sampling coherent 
and multi-thermal light fields, where the proposed model can be specified as follows. 
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3.A.I. Sampling a coherent light field 

For a coherent field, the photon-number distribution is given by 



2/i 



n,ph 



nl 



-\a\- 



11] 



for which n = 



(7 



(3) 



|a| 2 . Remembering that the distribution of dark-counts is also 



poissonian, according to Eq.s © and (fTUj) . we have: 



Q 



cl+dc 



1 



?™el+dc 
Qs, el+dc 



"T-el+dc 
- 3<5el+dc 



a 



(2) 

m,el+dc 
^el+dc 

1 



- 1 



a 



(2) 

m,el 



+ a {2) 

1 m,dc 



f m e i + m dc 



(m e i + m d 



m el+dc 



We thus obtain from Eq.s 



el+dc 

and (Til 



cr (3) 

u m, el+dc 
"^el+dc 



- 1 



a {3) 

w m,el 



a {3) 

u m,dc 



. m e i + m dc 



m e i + m dc 



x,out 
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x,out 



7- 



7 



3e 



7r 



1 + e 



(12) 
(13) 



at any mean value x out . 

The results in Eq.s ( JT2l) and ( fT3l) show that, if we measure a coherent field at different 
values of x on t by varying the overall detection efficiency r] through an attenuator [23] , we can 
obtain the values of 7, the detector gain, and e, the probability of cross-talk, by fitting the 
experimental values of F x>aat and S XtOUt . The 7- value is obtained in units of x per detected 
photon. Of course, no information can be obtained from Eq.s ([12]) and ( TT3|) on the amount 
of dark-counts as far as the statistics is poissonian for both light and dark-counts. Note that, 
as the relation between 7 and e is quadratic, we have in general two sets of possible values 
for each fit. The choice of one of the sets is made by comparing the value of 7 with the 
interspacing between the peaks of the pulse-height spectrum. 

3. A. 2. Sampling a multi-mode thermal light field 

The photon-number distribution of a field made of /i independent thermal modes each con- 



taining N t h/ fx mean photons 

p _ 

1 n,ph 
r(2) 



is given by 



{n + n-l)\ 



(14) 



n\(ji-l)\(N th /ti+l) tt ( t i/N th + l) n ' 

for which n = N th , ^ ] = N th (N th /Li + 1) and c^ 3 ) = N th (N th / > + 1) (2N th /fi + 1). In this 
case we have 



Q 



el+dc 



^el+dc 



""lel+dc 



(2) 

m, el+dc 



"^el+dc 



a 



(2) 
m,el 



(7 



(2) 
m,dc 



m e i + m dc 



(m e i + m dc ) 2 
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m c i(m e i//i + 1) + m dc - (m el + m dc ) 1 m, 



el 



(m e i + m dc ) 2 /i (m el + m dc ) 



Qs,el+dc 3Qel+dc 



2 

m el+dc 



1 


r<7 (3) 

"m,el+dc j 


a (2) 
q u m,el+dc 


+ 3 




2 

m el+dc 


"^el+dc 


TTTel+dc 




ff (3) , (3) _o (2) 
w m,el ^ °m,dc ou m,el 


°°m,dc i 


2 (mei + m dc ) 



(m e i + m dc ) 3 

m e i(m e i//i + l)(2m d //i + 1) + m dc - 3m e i(m e i//i + 1) - 3m dc + 2(m e] + m dc ) 



(m e i + m dc ) 3 

m 3 . 



/x 2 (m e i + m dc ) 3 



We now observe that the measured output mean value can be written as x out = 7(l+e)(m e i+ 
m dc ) = x p h + x dc , where only x p h undergoes attenuation during the experimental procedure 
while x dc remains constant. We thus rewrite Eq.s ([9]) and ( TTUj) as 

„ 1 A x dc \ 2 _ l + 3e 

/i V x out / 1 + e 

q 2 / x dc \ 3 _ 2 3 l+_3e / x dc \ 2 _ 2 l + 7e 

5 x ,out = —2(1 — = — I z ou t + -7-n — 1 - ^ z out + 7 -7- — • 16 

By measuring our multi-mode thermal field at different values of x out , and plotting F x ont 
and S^out w e can obtain an estimation of the parameters //, x dc , e, 7 from a self-consistent 
procedure and without need of an independent calibration. Again 7 is in units of x per 
detected photon and the relation between 7 and e is quadratic (see above). 

3.B. Model II: a numerical evaluation based on the photon-number resolving properties of 
SiPM 

The above mentioned self-consistent method is very powerful, but requires the acquisition of 
several histograms at varying 77, which could not always be easy to perform, or possible at 
all, in many practical applications: from this point of view, the possibility to analyze each 
spectrum independently looks complementary to the self-consistent approach. We performed 
this analysis with a two-step procedure: 

- we measured the areas of the spectrum peaks, thus obtaining an estimation of the 
number of counts per peak; 

- we fitted the obtained data points with a theoretical function, which takes into account 
the statistics of light, detection and all deviations of the detectors from ideality, such 
as DCR and cross-talk effects. 
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To evaluate the area of each peak, we performed a multi-peak fit of the spectrum histogram, 
modelling each peak with a Gauss-Hermite function |37j : 



GH = Ne~ w2/2 [l + 3 hH 3 (w) + 4 hH 4 (w)] , (17) 

where 

w = ^—^ (18) 
a 

and N is a normalization factor, x is the peak position and a is the variance of the gaussian 
function; H^(w) and H<±(w) are the third and the fourth normalized hermite polynomials 
and their contribution gives the asymmetry of the peak shape, whose entity is regulated by 
the pre-factors 3 h and 4 h, with values in the range [—1,1]. The global fit function of the 
spectrum is a sum of as many Gauss-Hermite function as the number of resolved peaks. 
The choice of the GH-function in Eq. (fT7j) allows us to calculate the area A n of the n-th 
peak in a very straightfoward way, simply by the relation 

A n = N n a n (V2^+ 4 h n ). (19) 

The error OA n on the obtained value is calculated by propagating the errors on the fit 
parameters. 

This analysis is also useful in order to calculate the system gain 7: in fact, from the fitted 
values of the peak positions x n we can calculate the peak-to-peak distance A for all the 
resolved peaks: 

Ara,n+1 = X~n+1 ~ %n- (20) 

The error a& n n+l associated to this value is once again obtained by propagating the fit errors 
of the peak position values; furthermore, to estimate 7 a weighted average on all the peak- 
to-peak values obtained from the analyzed histogram is performed. 

As for the theoretical function, the effect of detection, DCR and amplification is modelled 
as described in previous sections (see Eq.s ©-flH]) and Eq. (JS])). 

The effect of cross-talk, is described by using a bernoullian process, in a way analogue to what 
has been done with function Ck, m ( e ) of Eq. ([6]). However, as cross-talk process is intrinsically 
a cascade phenomenon, its contribution has been calculated by adding higher order effects: 

k m n 

Pk,cross ^ ] ^ ] ^ ] Pk—m—n—j,el+dc^m,k—m—n—j(^^n,rn(^}^j,n(^)i (^1) 
m=0 n=0 j=0 

where terms like Bj^ n (e) stand for the bernoullian distribution 

BjM = I k'(l-e)^. 
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Such a higher order expansion is not trivial to be achieved by the self-consistent approach of 
Model I, in which an explicit analytic expression of P x ,out is needed in order to calculate its 
momenta. Here, as all the elements of interest (m e i, m^c, e, the number of modes fi) will be 
obtained as fit parameters, this is not necessary and therefore P Xj0 ut can be just numerically 
calculated as the fitting function. 

The major limitation of this approach is obvious: as all the information on the statistics 
of the system is obtained from the peak areas, this method can only be applied to peak- 
resolving histograms with a number of peaks greater than the number of free parameters of 
the fitting function, which, in the present analysis, can rise up to five. 

4. Experimental results 

4- A. Coherent light 

First of all, we measured the coherent light emerging from the laser. We measured the values 
of the output, x, at 20000 subsequent laser shots for 15 series, each one with a different mean 
value x ou t, set by rotating a polarizer (P in Fig. [3]) in front of the collection fiber. We then 
acquired a series of data in the absence of light and set the zero at the mean value of the 
main peak of the resulting histogram, in which the presence of dark counts and cross talk 
emerges as a much lower separated peak. Following Model I we evaluated the experimental 
values of F Xt0ut and S^out that are plotted in Fig. H] as a function of x out . We then fitted 
the data to straight lines and obtained (81.1 ± 0.2) ch, for -F x out , and (6971 ± 57) ch 2 , for 
S x ,out- These values were used to evaluate 7 and e from Eqs. (fl2]) and (|T3|) . We obtained 
7 = (75.4 ± 1.3) ch and e = (0.039 ± 0.009). The x- values were then divided by 7 and 
re-binned in unitary bins j23J[3B] to obtain the Pfc jC ross distribution of the actual light signal 
amplitude measured in the presence of dark-counts and cross-talk. Note that due to the lin- 
earity of the detector the mean value of the output can be directly obtained as k CTOSS = x out /7, 
independent of the shape of the distribution. 

In Fig. |5] we plot as bars six different Pfc,cross distributions at different mean values. Super- 
imposed to the experimental values we plot two theoretical curves, one is a poissonian (see 
Eq. ( ITT]) ) having mean value k cross (white circles), while the other (full circles) is evaluated by 
including the cross-talk effect. We evaluate Pfc.cross from Eq. ((7j) in the case in which P mie i + d c 
is poissonian and get: 



\™. = e-^-Hl - r F, (l, I - * 1 - f ; J 1 -'?"** ) E^M , (23) 



where m c i +dc = k CTOSS /(l + e) = x out /(7(l + e)). The theoretical distributions (full circles in 
Fig. [5]) are evaluated by using the measured values of x out , 7 and e. 

Note that measuring a coherent light with this method enables the simultaneous character- 
ization of the detector gain and of the contribution of cross-talk. The comparison between 
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the data and the theoretical functions can be estimated through the evaluation of the fidelity 

m 

/ — X/ V -Pfc,exp-Pfc,theo • (24) 
k=0 

On the other hand, using the method of analysis of Model II, we could study each one of 
the acquired histograms separately. As mentioned above, in the case of coherent light, we 
have a theoretical fitting function with a total of three free parameters: expectation value 
of light and DCR contribution m e i + d c , the probability e for an avalanche to trigger a second 
one and a global normalization factor (up to three "iterations"): in this case we are limited 
to spectra with at least 4 resolved peaks. 

In Fig. [6] we show the results of this analysis for the multi-peak fit and for the fit of the 
statistics of the avalanches using \ 2 as an indicator of the goodness of the fit result. The 
obtained e values are compatible with what we obtained with Model I. Also the gain values, 
evaluated as the peak-to-peak distance, show a good agreement with the 7 values given by 
Model I: 7 = (78.28 ± 0.26) ch for spectrum in left panel and 7 = (72.61 ± 0.19) ch for 
spectrum in right panel. 

4-B. Multi-mode pseudo-thermal light 

In order to obtain information on the contribution of dark-counts we have to measure a 
different light statistics, whose shape is modified by the convolution with the poissonian 
distribution for dark-counts. We thus produced a pseudo-thermal light field by selecting 
with a small aperture (~ 150 fim diameter) a region much smaller than the coherence area 
of the speckle patterns produced by the rotating diffuser. We follow the same procedure 
described for coherent light by measuring the values of the output, x, at 50000 subsequent 
laser shots and at 10 different mean values, obtained by means of a variable neutral-density 
filter (ND in Fig. [3]). In Fig. [7] we plot the experimental values of F x out and S x>out as a function 
of x out . Along with the experimental data we plot the fitting curves evaluated according to 
Eqs. ( fT5l) and ( fl6l) . To describe the fitting procedure in detail we rewrite Eqs. ( fl5l) and ( JT6l) 
as the 

F x , out = (l-^-Yxout + B (25) 

V ^out / 

5,, out = A(l-^) 3 xl nt + 3B(l-^-) 2 x out + C , (26) 

\ -^out / \ ^out / 

where we have set /x = 1. First of all we fitted the data to F x>out and obtained the values of 
x dc = (5.82028 ± 1.34015) and B = (87.805 ± 2.09009) ch. Then we fitted the data for S x>out 
by substituting the obtained values of x^c ch and B to obtain A = (2.34754 ± 0.091576) 
and C = (8531.48 ± 419.571) ch 2 . These values are then used to evaluate 7 and e from 
Eqs. ([I5D and (d6]). We obtained 7 = (74.2785 ± 18.6017) ch and e = (0.100174 ± 0.166894). 
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The x-values were then divided by 7 and re-binned in unitary bins [23] to obtain the -P^cross 
distribution of the actual light signal amplitude measured in the presence of dark-counts 
and cross-talk. In Fig. [8] we plot, as bars, six different Pk, C mss distributions at different mean 
values. Superimposed to the experimental values we plot two theoretical distributions: the 
first one (open circles) is evaluated by including the contribution of dark-count that modifies 
the statistics of a single-mode thermal distribution (see Eq. (|14p ) into P m)C i+dc according to 
Eq. f[2"Tj) . which, in the present case, yields: 



Pm,c\+Ac ^ ] Pk,dcPm—k,c\ 

k=0 



1 + JL^-d + TiYu 



a 

— m, 1 — m — ll } m<jc 1 + 



iP) 



(ji - 1)! V m e J \ 11 J 

where U (a, b, z) is the confluent hypergeometric function. The parameters are evaluated as 
= x&d (7(1 + €.)) and m e \ = (x out — Xdc) / (7(1 + The second distribution (full circles) 
is evaluated from Eq. (J7J) to take into account the cross-talk. Unfortunately, the calculation 
does not yield an easy analytical result, and hence we evaluate it numerically. The values of 
the fidelity for the data in Fig. |H] improve when we take into account both dark-counts and 
cross-talk. 

Turning now to the other approach, we note how the number of fit parameters in this case 
is enhanced: we now have the expectation value m e i of avalanches generated by detection, 
the expectation value md c of DCR contribution, the number of modes /1, the probability e of 
triggering a cross-talk event (up to three iterations) and again a global normalization factor, 
for a total of 5 fit parameters: obviously, this puts a severe limit on the applicability of this 
method, needing at least 6 resolved peaks. 

As it can be noted from the fit results in Fig. [UJ once again the results obtained by using 
Model II are compatible within errors with what we found by applying Model I. However, 
whether the global fits present a very low x 2 value for degree of freedom, the obtained fit 
parameters present high uncertainties, probably indicating the presence of very high off- 
diagonal elements in the minimization matrix and suggesting a strong correlation between 
the various parameters. This problem can be avoided by fixing some of the fit parameters 
(such as DCR or cross-talk), by retrieving their value from an accurate direct measurement 
as that explained in Section [2j 

5. Discussion 

We can compare the results of the two analysis methods on the same data-sets taken at 
different intensities. In Fig. [10] (a) we plot as full circles the values of e obtained by applying 
Model II to coherent light along with their weighted average (full line). As a comparison, 
the value of e obtained by applying Model I is plotted as dashed line. In Fig. [10] (b) we plot 
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as full circles the values of mean photon numbers evaluated for the same data as in panel (a) 
by applying Model II. White circles represent the values of mean photon numbers obtained 
by applying Model I. As we can see, the values are compatible within errors. 

The same comparison for the measurements on thermal light is shown in Fig. [TTJ Here 
we can see that the agreement is better for the mean values of detected photons (panel (a)) 
and for the DCR (panel (b)), while the estimated values of e from the two Models definitely 
disagree. This can be due to the different approximations adopted by the two Models (first 
order vs third order) that become relevant when measuring thermal light instead of coherent 
light. 

In Table |2] we summarize the results of the two Models. We demonstrated that both the 
Models work in a self-consistent way, even if they have two definitely different approaches. 
Model I does not need peak resolving capability, but requires the acquisition of several 
histograms at varying r\. Once determined the parameters e and DCR, all the data-sets in a 
series can be analyzed, independent of the number of distinguishable peaks in the pulse-height 
spectrum. Model II works analyzing each histogram independently, but, as GM- avalanches 
distribution is obtained with a fit of the data, it requires at least a number of resolved peaks 
greater than the number of free parameters. The fact that the two Models give very similar 
results for mean photon numbers is particularly important as in most applications this is the 
only important parameter. Merging the two Models we can devise an optimal strategy based 
on a self-consistent calibration performed by measuring a known light and analyzing the 
data with Model I: once known e and DCR, the determination of the mean photon number 
is independent of the specific statistics of light. Hence the information on the mean photon 
number can be obtained from each single measurement, even when the fitting procedure of 
Model II cannot be applied. 

6. Conclusions 

We have developed a model to describe the operation of a multi-pixel detector for which 
dark-count and cross-talk effects are non negligible and we have implemented two different 
procedures for recovering the values of dark-counts and cross-talk, both implemented on 
measurements performed on the same light under characterization. We demonstrated that 
Model I is self-consistent and does not need peak resolving capability; on the other hand, it 
requires the acquisition of several histograms at varying rj. Model II works by analyzing each 
histogram independently, but seems to be a less robust method in order to perform direct 
measurements of the detector parameters and it is better performing if supported with direct 
measurement of some relevant parameters. The results obtained by both methods show that 
the light statistics can be reliably reconstructed in the case of coherent and thermal light. 
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Fig. 1. Left panel: Output of a SiPM Hamamatsu MPPC S10362-11-100C detecting a weak 
light field as displayed by a 100 MHZ oscilloscope. Right panel: histogram of the correspond- 
ing spectrum. 

Fig. 2. Threshold scan of the SiPM performed at room temperature with no impinging light. 
Fig. 3. (Color online) Experimental setup. Nd:YLF: laser source, P: polaroid, ND: variable 
neutral density filter, SiPM: detector. The components in the dashed boxes are inserted to 
produce the pseudo-thermal field. 

Fig. 4. (Color online) Plot of i^ j0u t an d S^out as a function of x out for coherent light. 
Fig. 5. (Color online) Experimental i\ C ross distributions at different mean values (bars) and 
theoretical curves evaluated according to Model I: poissonian (white circles), poissonian mod- 
ified by cross-talk effect (full circles). The corresponding fidelity values of the reconstruction 
are also shown. 

Fig. 6. (Color online) Experimental results for Model II applied on two of the histograms 
acquired with coherent light. Upper row: result of the multi-peak fit procedure; lower row: 
fitted theoretical distributions. The corresponding fidelity values of the reconstruction are 
also shown. 

Fig. 7. (Color online) Plot of F x>out and S X)OUt as a function of x out for pseudo-thermal light. 
Fig. 8. (Color online) Experimental -P^cross distributions at different mean values (bars) 
and theoretical distributions evaluated according to Model I: thermal modified by dark 
count distribution (white circles), thermal modified by dark counts and cross-talk effect (full 
circles). The corresponding fidelity values of the reconstruction are also shown. 
Fig. 9. (Color online) Experimental results for Model II applied on two of the histograms 
acquired with thermal light. Upper row: result of the multi-peak fit procedure; lower row: 
fitted theoretical function. The corresponding fidelity values of the reconstruction are also 
shown. 

Fig. 10. (Color online) Left panel: values of e obtained by applying Model II to coherent 
light (full circles) and their weighted average (full line). Dashed line: value of e obtained by 
Model I. Right panel: values of mean photon numbers evaluated by applying Model II (full 
circles) and by applying Model I (white circles). 

Fig. 11. (Color online) Left panel: values of e obtained by applying Model II to thermal 
light (full circles) and their weighted average (full line). Dashed line: value of e obtained 
by Model I. Central panel: values of DCR evaluated for the same data as in the left panel 
by applying Model II (full circles) and their weighted average (full line). Dashed line: value 
of DCR obtained by Model I. Right panel: values of mean photon numbers evaluated by 
applying Model II (full circles) and by applying Model I (white circles). 
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Fig. 1. Left panel: Output of a SiPM Hamamatsu MPPC S10362-11-100C 
detecting a weak light field as displayed by a 100 MHZ oscilloscope. Right 
panel: histogram of the corresponding spectrum. 
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Fig. 2. DCR threshold scan of the SiPM performed at room temperature with 
no impinging light. 
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Fig. 3. (Color online) Experimental setup. Nd:YLF: laser source, P: polaroid, 
ND: variable neutral density filter, SiPM: detector. The components in the 
dashed boxes are inserted to produce the pseudo-thermal field. 
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Fig. 4. (Color online) Plot of F X)OUt and S XtOVLt as a function of x out for coherent 
light. 
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Fig. 5. (Color online) Experimental -Pfc iCroS s distributions at different mean 
values (bars) and theoretical curves evaluated according to Model I: poisso- 
nian (white circles), poissonian modified by cross-talk effect (full circles). The 
corresponding fidelity values of the reconstruction are also shown. 
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Fig. 6. (Color online) Experimental results for Model II applied on two of the 
histograms acquired with coherent light. Upper row: result of the multi-peak 
fit procedure; lower row: fitted theoretical distributions. The corresponding 
fidelity values of the reconstruction are also shown. 
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Fig. 7. (Color online) Plot of F X)OUt and S XtOUt as a function of af ou t for pseudo- 
thermal light. 
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Fig. 8. (Color online) Experimental -P^cross distributions at different mean 
values (bars) and theoretical distributions evaluated according to Model I: 
thermal modified by dark count distribution (white circles), thermal modified 
by dark counts and cross-talk effect (full circles). The corresponding fidelity 
values of the reconstruction are also shown. 
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Fig. 9. (Color online) Experimental results for Model II applied on two of the 
histograms acquired with thermal light. Upper row: result of the multi-peak 
fit procedure; lower row: fitted theoretical function. The corresponding fidelity 
values of the reconstruction are also shown. 
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Fig. 10. (Color online) Left panel: values of e obtained by applying Model II to 
coherent light (full circles) and their weighted average (full line). Dashed line: 
value of e obtained by Model I. Right panel: values of mean photon numbers 
evaluated by applying Model II (full circles) and by applying Model I (white 
circles) . 
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Fig. 11. (Color online) Left panel: values of e obtained by applying Model II to 
thermal light (full circles) and their weighted average (full line). Dashed line: 
value of e obtained by Model I. Central panel: values of DCR evaluated for 
the same data as in the left panel by applying Model II (full circles) and their 
weighted average (full line). Dashed line: value of DCR obtained by Model I. 
Right panel: values of mean photon numbers evaluated by applying Model II 
(full circles) and by applying Model I (white circles). 
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Hamamatsu MPPC S10362-11-100C 



Number of Diodes: 100 

Area: 1 mm x 1 mm 

Diode dimension: 100 /xm x 100 /xm 

Breakdown Voltage: 69.23 V 

dark-count Rate: 540 kHz at 70 V 

Optical Crosstalk: 25 % at 70 V 

Gain: 3.3 • 10 6 at 70 V 

PDE (green): 15 % at 70 V 



Table 1. Main characteristics of the SiPM (Hamamatsu, model MPPC S10362- 
11-100C). The data refer to room temperature. 
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Poissonian 


Pseudo-thermal 




Model I 


Model II 


Model I 


Model II 




DCR 






0.071 ± 0.027 


0.094 ± 0.035 




e 


0.038 ± 0.008 


0.027 ± 0.005 


0.091 ± 0.138 


0.035 ± 0.004 



Table 2. Comparison between the global DCR and cross-talk values obtained 
with Model I and weighted average of the values obtained with Model II. 
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